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ABSTRACT 


f  The  dynamical  response  of  a  three-dimensional  hydroelastic 
model  of  the  cochlea  is  studied  for  a  pure  tone  forcing.  The 
basilar  membrane  is  modeled  as  an  inhomogenous,  orthotropic 
elastic  plate  and  the  fluid  is  assumed  to  be  Newtonian.  The 
resulting  mathematical  problem  is  reduced  using  viscous  boundary 
layer  theory  and  slender  body  approximations.  This  leads  to  a 
nonlinear  eigenvalue  problem  in  the  transverse  cross-section. 

The  solutions  for  the  case  of  a  rectangular  and  semi-circular 
cross-section  are  computed  and  comparison  is  made  with 
experiment.  The  role  of  the  place  principle  in  determining  the 
difference  limen  is  presented  and  it  is  shown  how  the  theory 
agrees  with  the  experimental  measurements. 
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INTRODUCTION 

The  theoretical  and  experimental  basis  for  our  under¬ 
standing  of  the  mechanical  resolution  of  sound  by  the  ear, 
principally  as  takes  place  in  the  cochlea,  is  currently 
undergoing  an  extensive  reexamination.  To  understand  why,  it  is 
necessary  to  outline  how  the  basic  system  functions.  After  the 
signal  is  transmitted  through  the  outer  and  middle  ear  it  is 
received  by  the  cochlea  through  the  movement  of  the  stapes, 
which  is  the  last  bone  in  the  ossicular  chain.  This  motion 
causes  the  fluid  inside  the  cochlea  to  move  which,  in  turn, 
stimulates  waves  to  propagate  down  the  basilar  membrane 
(Fig.  1).  This  membrane  widens,  as  well  as  decreases  in 
thickness,  so  that  the  amplitude  of  the  wave  begins  to  grow  and, 
at  the  same  time,  the  wavelength  decreases.  However,  because  of 
damping,  the  wave  eventually  reaches  a  point  on  the  membrane 
after  which  it  quickly  decays  to  zero.  The  location  and 
absolute  level  of  this  maximum  amplitude  depends  on  the 
frequency;  the  higher  the  frequency  the  closer  the  location  is 
to  the  stapes.  The  result  of  this  dependence  is  a  frequency  map 
on  the  basilar  membrane.  It  is  thought  that  this  leads  to  the 
mechanical  resolution  of  the  signal  since  there  are  sensory 
cilia  setting  on  the  membrane  which  respond  to  its  deflection 
and  transmit  neural  signals  accordingly. 

The  problem  with  this  form  of  mechanical  resolution  is  that 
it  does  not  seem  to  be  sharp  enough  to  explain  the  extremely 
fine  tuning  that  is  observed  in  the  neural  signals  leaving  the 
cochlea.  This  conclusion  comes  from  experimental  (Kiang  et  al. 


1965;  Geisler  et  al,  1974;  Russell  and  Sellick,  1978)  as  well  as 
theoretical  (Zwislocki,  1965;  Steele  and  Taber,  1979;  Holmes, 
1982)  analysis  of  the  system.  It  is  for  this  reason  that  there 
has  been  a  search  in  the  last  few  years  for  a  "second  filter"  to 
explain  the  fine  tuning,  but  as  of  yet  they  have  been 
unsuccessful  (Frommer,  1979;  Zwislocki,  1980;  Lighthill,  1981; 
Khanna,  1983). 

The  problem  is  made  more  difficult  by  the  current  state  of 
the  experimental  evidence  that  is  available.  Until  recently  some 
of  the  best  results  were  obtained  using  Mossbauer  measurements  of 
the  basilar  membrane's  motion  (Johnstone  and  Boyle,  1967;  Rhode, 
1971).  The  procedure  involves  opening  the  cochlea,  placing  a 
radioactive  source  on  the  basilar  membrane,  and  then  measuring 
the  response  over  several  hours.  Unfortunately,  it  appears  that 
the  radioactivity  has  a  detrimental  effect  on  the  viability  of 
the  cells  in  a  relatively  short  time  (Kliauga  and  Khanna,  1983), 
and  so,  the  results  of  these  experiments  are  suspect.  Moreover, 
the  cochlea  has  an  intricate  defense  mechanism  and  is, 
apparently,  very  sensitive  to  outside  interference  (Khanna  and 
Leonard,  1981  and  1983).  This  means  that  the  available 
experimental  techniques  used  to  date,  including  the  laser 
interferometric  measurements  of  Khanna  and  Leonard  (1982),  are 
probably  only  able  to  give  a  gross  description  of  the  normal 
response  and  are  not  capable  of  measuring  the  very  fine  tuning, 
if  present,  that  is  sought. 

It  is  with  this  rather  unsettling  state  of  affairs  that  we 
consider  the  three-dimensional  hydroelastic  model  for  the 


cochlea  that  is  outlined  below.  Even  though  a  considerable 
amount  o£  work  has  been  done  on  it  (Steele,  1976;  Chadwick  and 
Cole,  1979;  Steele  and  Taber,  1979;  Holmes,  1980;  Holmes,  1982) 
it  has  never  been  completely  solved  and  probably  never  will 
except  by  tour  de  force  numerical  methods.  It  is  therefore 
necessary  to  use  certain  approximations,  such  as  boundary  layer 
theory,  to  reduce  the  problem  to  a  more  tractable  size.  Steele 
(1974,  1976,  1980)  has  made  significant  contributions  to  the 
resolution  of  this  problem  and  they  have  been  the  basis  for  much 
of  our  own  work.  His  approach  involves  a  WKB  type  of 
approximation  of  the  solution,  in  which  he  specifies  the  shape 
of  the  transverse  displacement  of  the  basilar  membrane.  With 
this,  energy  methods  and  time  averaged  Lagrangians  are  used  to 
find  the  eikonal  and  transport  equations.  Our  approach,  on  the 
other  hand,  uses  perturbation  expansions  of  the  original  dif¬ 
ferential  equations  and  systematic  viscous  boundary  layer 
theory.  Chadwick  (1980)  and  Holmes  (1982)  have  studied  the 
fully  three-dimensional  problem  in  this  way,  but  in  doing  so 
other  assumptions  such  as  high  frequencies  or  a  small  mass 
density  for  the  basilar  membrane  are  imposed.  We  intend  in  this 
article  to  remove  these  restrictions  and  obtain  the  leading 
order  approximation  to  the  full  hydroelastic  problem  for  a  pure 
tone  forcing.  As  with  any  mathematical  formulation  of  a  complex 
physical  problem,  various  assumptions  are  made  concerning  the 
structure  and  composition  of  the  system,  but  they  include  most 
of  these  earlier  models. 
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FORMULATION  OF  MODEL 

We  will  assume  the  geometrical  structure  of  the  cochlea  to 
be  such  that  it  consists  of  an  unrolled  tapered  tube  containing 
two  fluid  chambers.  The  partition  separating  these  chambers 
consists  of  a  rigid  portion  (the  bony  shelf),  a  flexible  portion 
(the  basilar  membrane),  and  an  aperture  at  the  apical  end  (the 
helicotrema ) .  The  variability  of  the  geometry  is  arbitrary 
except  that  we  will  assume  the  cochlear  wall  to  be  symmetric 
through  the  partition,  in  this  case,  through  the  x,y  plane.  The 
cochlear  wall,  along  with  the  shelf,  are  rigid  except  at  the 
basal  end  where  there  are  two  openings  (the  oval  and  round 
windows)  that  are  covered  by  flexible  membranes.  The  footplate 
of  the  stapes  is  attached  to  the  cochlea  at  the  oval  window. 

Although  the  assumptions  made  on  the  geometry  are  not 
unreasonable,  they  do  impose  limitations  on  the  theory.  The 
spiraling  of  the  cochlea  is  neglected  as  it  is  believed  to 
effect  the  response  at  only  the  lowest  frequencies  and  is 
primarily  for  "packaging  purposes  "  Also,  the  two  chambers  are 
not  exactly  symmetric  and  there  is  a  third  chamber,  the  scala 
media,  which  accounts  for  about  8%  of  the  total  volume.  The 
possible  effect  of  this  third  chamber  has  been  studied  by 
Adjemian  (1981)  and,  although  the  results  are  incomplete,  it 
does  not  appear  to  have  a  significant  effect  on  the  mechanical 
response  of  the  system.  It  is,  however,  important  for  the 
neural  response  because  of  the  different  cations  that  these 


chambers  contain 


OVAL  WINDOW 


FIG.  1.  The  geometry  and  notation  used  for  the  hydroelastic  model  of  the  cochlea. 
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The  fluid,  called  perilymph,  that  fills  the  two  chambers  is 
similar  in  terms  of  its  mechanical  response  to  water  (Money  et 
al,  1966;  Rauch  and  Rauch,  1974).  Consequently,  it  is 
considered  to  be  Newtonian.  Compressibility  is  omitted,  even 
though  the  frequencies  in  the  audible  range  can  extend  up  to 
15,000  Hz,  as  the  overall  length  of  the  cochlea  is  only  about 
3.5  cm  and  the  pressures  at  the  stapes  are  relatively  small 
(Lighthill,  1981).  The  fluid  in  the  scala  media,  called  endo- 
lymph,  contains  a  relatively  large  concentration  of  proteins 
(Schuknecht,  1970;  Rauch  and  Rauch,  1974).  These  macromolecules 
could  very  well  have  a  significant  effect  on  the  response  but 
are  neglected  in  our  study.  Also,  the  convective  terms  in  the 
Navier-Stokes  equations  are  not  included  in  our  analysis.  The 
reason  is  that,  after  scaling  the  problem,  one  finds  that  these 
nonlinear  terms,  relative  to  those  obtained,  are  on  the  order  of 
the  ratio  of  the  stapes  amplitude  to  the  jross -sectional  width 
of  the  cochlea.  For  normal  sound  pressures  this  ratio  is 
around  10~4  and  so  they  do  not  contribute  to  the  first  term. 

As  for  the  basilar  membrane,  it  is  a  composite  material 
made  up  of  an  orthogonal  weave  of  fibers  and  interspersed  with 
ground  substance  (Iurato,  1967).  From  the  studies  of  its 
elastic  properties  it  appears  to  respond  mechanically  as  an 
elastic  plate  (von  Bekesy,  1960;  Novoselova,  1975).  Therefore, 
because  of  its  fiber  network  and  its  variable  thickness  we  will 
model  it  as  an  inhomogenous  orthotropic  elastic  plate.  It  is 
not  necessary  to  specify  its  exact  location  and  configuration  in 
the  partition  at  this  time  as  the  analysis  to  follow  applies  to 
essentially  any  shape  basilar  membrane. 


The  remaining  approximation  concerns  the  relatively  high 
frequencies  in  the  auditory  range  (about  25  Hz  or  higher  for  the 
human  cochlea).  The  Reynolds  number,  which  is  based  on  the 
width  of  the  cross-section,  is  relatively  large,  and  so, 
boundary  layer  theory  can  be  used  for  the  fluid.  The  error  in 
this  approximation  is  related  to  the  curvature  of  the  wall.  In 
the  examples  to  be  considered  later  there  are  corners,  but  since 
they  are  isolated  they  do  not  interfere  with  the  validity  of  the 
boundary  layer  approximation  (Holmes,  1981). 

The  pressure  of  the  fluid  is  represented  by  p(x,y,z,t)  and 
the  vertical  displacement  of  the  basilar  membrane  is  n(x,y,t). 

In  what  follows  both  the  independent  and  dependent  variables  are 
dimensionless,  but,  where  necessary,  an  asterisk  is  used  to 
indicate  their  dimensional  analogue.  So,  for  example,  the 
nondimens ional  spatial  coordinates  x,y,z  are  related  to  x*,y*,z* 
as  follows 


x*  =  Lx  ,  y*  =  By  ,  z*  =  Bz  , 


where  B  and  L  are  the  characteristic  width  and  length  of  the 
basilar  memembrane  respectively.  As  we  are  interested  in  the 
response  to  a  pure  tone  let 


and 


p(x,y,z,t)  =  p(x,y,z)eifc  , 


n(x,y,t)  *  ;(x,y)eifc 


wh  re  t  -  .t*  and  Q  is  the  given  driving  frequency.  Also  note 
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p*(x*,y*,z*,t*)  =  pAs*Lft2p(x,y,z,t) 

and 

n*(x*,y*,t*)  =  ^pAs*L^  n(x,y,t)  , 

where  As*  is  the  amplitude  of  the  stapes,  p  is  the  fluid  den¬ 
sity,  and  uc  is  the  mass  per  unit  area  of  the  basilar  membrane 
at  its  center  (x=l/2).  The  pressure,  because  of  the  symmetry, 
is  an  odd  function  of  z,  and  so,  it  is  only  necessary  to  con¬ 
sider  the  fluid  motion  in  the  upper  chamber  (f). 

From  the  above  discussion  we  therefore  have  the  following 
equations  of  motion: 


i  )  for  the  fluid 

U23x2  +  3y2  +  3z2)p  =  0  ,  (la) 

ii)  for  the  basilar  membrane 

3y2(D23y2C)  +  e2[ 3y2( V1D23X2C )  +  3X2 ( vl°2 8y2 t )  +  4 3X SyfD^ 3X 3y 4)  J 

+  e43x2(D13x2;)  -  ui2h  c  =  -2u>2p(x,y,0+) .  (lb) 


The  parameters  are 

R  ,  i) 

e  =  £  and  w  =  _  , 

L  nc 

where  ac  represents  a  characteristic  frequency  of  a  beam  in  the 
cross-section  and  is  given  by 


a  2  =  ^2c* 
c 


Here  D2c*  is  the  bending  rigidity  in  the  y-direction  at  the  cen¬ 


ter  of  the  basilar  membrane.  The  functions  ,  d2,  and  are 


the  respective  bending  and  twisting  rigidities  normalized  by 
D2c*»  and  \>i  is  the  Poisson  ratio  characterizing  the  decrease 
in  the  y-direction  during  tension  applied  in  the  x-direction 
( Lekhnitskii ,  1968).  Also,  h(x,y)  represents  the  thickness 
of  the  basilar  membrane  after  normalization  by  hc,  which  is  the 
thickness  at  the  center. 

The  kinematic  boundary  condition  for  the  pressure 
(Appendix  A),  when  0<x<l,  is 

I  -a c  on  the  basilar  membrane 

(2) 

0  on  the  rigid  wall  , 

where  n  represents  the  unit  outward  normal.  The  6  term  repre¬ 
sents  the  viscous  boundary  layer  correction  to  the  inviscid 
boundary  condition,  where  the  constant  0  is  given  as 


a  =  Pi 
yc 

and  v  is  the  kinematic  viscosity  of  the  fluid. 

It  is  not  necessary  at  the  moment  to  specify  the  boundary 
conditions  for  the  plate,  but  we  will  assume  that  it  is  either 
clamped  or  simply  supported.  Also,  the  boundary  conditions  at 
the  basal  and  apical  ends,  where  x  =  0,1,  are  discussed  later. 
It  therefore  remains  to  solve  the  above  coupled  system  of 
partial  differential  equations  to  be  able  to  describe  the 
mechanical  response  to  a  pure  tone. 
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TWO-VARIABLE  EXPANSION 

For  most  mammalian  cochlae  the  aspect  ratio  e  is  very 
small.  For  example,  for  humans  it  is  on  the  order  of  10“  .  We 
can  take  advantage  of  this  through  the  use  of  a  slender  body 
approximation.  For  the  case  at  hand,  one  finds  that  the 
appropriate  expansions  are 

p  -  ee"^k^x»e^  [pQ(x,y,z )+epi+- • • ]  ,  ( 3a 


C  -  ee 


-ik(x,e)  [ 50(x,y)+eci+. . • ]  , 


where 


k(x, e  )  =  e ( x  y 


Introducing  these  into  (1),  the  0(1)  problem  is  found  to  be 

( 3y2+3z2)P0  =  ©x2P0  * 

and 

dy2 ( D23y2C o )  ”  Qx2[3y2<  vlD2C0>  +  vlD23y2«0+43y(Dk3y^0>  J 

+  (Qx4D1  -  u>2h)e0  =  -2o)2p0(x,y,0+) 

As  for  the  kinematic  boundary  condition  (2),  note  that  the 
boundary  of  the  upper  chamber  can  be  parameterized  as 
y  =  y (s , x ) ,  z  »  z(s,x)  where  0<s<l  represents  the 
cross-sectional  variation  and  0<x<l  the  longitudinal.  Thus 
the  normal  derivative  in  (2)  expands  as  follows 


3n  *  3nT  +  £2b<x's>I3x  -  \  b<x,s>3nT]  +  0( e4)  ' 


where  n<r  is  the  unit  outward  normal  in  the  transverse 
cross-plane  (x  fixed)  and 


b(x,s ) 


With  this  and  (2)  we  have  that 


_  i  -a^Q  on  BM 
(  3 n  ~03n  2)Po  =  , 

T  T  (  0  on  rigid  wall 


(5) 


To  be  definite  we  will  now  assume  that  the  plate  is  simply 
supported,  and  so, 

Co  =  3y2C 0  =  0  for  y  =  G+(x),G_(x)  .  (6) 

The  boundary  y  =  G+(x)  represents  the  portion  attached  to  the 
spiral  lamina  and  y  =  G_(x)  is  the  portion  attached  to  the 
spiral  ligament.  It  should  be  emphasized  that  clamped 
conditions,  or  a  combination  of  simply  supported  and  clamped 
boundary  conditions  could  also  be  used. 

The  0(1)  problem  is  now  essentially  complete.  From  it,  it 
is  only  possible  to  determine  uniquely  ©x2  along  with  the  y,z- 
dependence  of  Co  an<3  PO*  it  does  not  resolve  the  x-dependence 
of  these  variables  and  to  do  this  one  needs  to  consider  the  0(  e) 
problem.  The  details  are  outlined  in  Appendix  B,  and  they 
basically  involve  integrating  this  higher  order  problem  over  the 
cross  section.  To  simplify  the  expression  we  will  consider  the 
case  of  when  the  bending  and  twisting  rigidities  are  independent 
of  y,  and  are  only  functions  of  x.  The  result  is  that  pq  and 
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Cq  have  to  be  such  that 


where 


D3(x)  =  v^D2  +  2D){  . 


To  complete  the  problem  we  now  need  to  consider  the 
boundary  conditions  at  x  =  0  and  at  x  =  1.  At  the  basal  end  we 
have,  using  a  boundary  layer  approximation  as  in  (2),  that 

_  (  "n  on  oval  window 

Ox+£63x2)p  =  (8) 

(  0  on  rigid  portion, 

where  "n(y,z)  represents  the  known  displacement  of  the  stapes 
footplate.  To  use  this  with  the  two-variable  expansion  (3)  it 
is  necessary  to  introduce  an  edge  layer,  with  coordinate 
x  =  x/e,  and  match  the  expansions  in  this  region  with  (3).  One 
finds  that  for  the  two  expansions  to  match  to  the  first  order, 
mass  must  be  conserved,  and  so,  at  x  =  0  we  require 


and 


ex(l-ieex)  J  fpQ  dydz  =  i / fH  dydz  , 
*  rw 


k  (  0, e  )  =  0  . 


(9) 


(10) 


This  omits  the  contribution  of  the  basilar  membrane  which, 
because  of  its  tapered  geometry,  is  assumed  to  be  negligible 
compared  to  the  movement  of  the  fluid  in  this  region.  Of 
special  interest  is  the  case  of  when  "n  =1,  so  that  at  x  =  0 
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we  require 


/ / p«  dydz  =  w  ,  (11) 

JiP°  y  ex(i-i eex) 

where  A*  represents  the  (normalized)  area  of  the  round  window. 

It  would  appear  from  (7)  that  there  is  only  one  constant  of 
integration.  However,  there  are  at  least  two  in  the  general 
solution  since  there  are  two  values  of  k(x,e)  that  are  obtained 
from  the  square  root  of  0X2.  In  fact,  the  general  solution 
consists  of  the  superposition  of  the  eigenfunctions  obtained 
from  (4).  It  is  therefore  necessary  to  introduce  an  edge  layer 
at  x  «  1.  However,  by  assuming  the  continuity  of  the  pressure, 
and  since  it  is  an  odd  function  of  z,  it  follows  that  at  the 
apical  end 


PodrYfZ)  ■  0 


With  this  last  boundary  condition  the  problem  for  the  first 
term  approximation  in  (3)  is  essentially  complete.  It  is 
centered  around  the  unusual  eigenvalue  problem  given  in  (4), 
is  the  eigenvalue  and  pg,  eg  are  the  corresponding  eigen¬ 
functions.  The  problem  is  a  nonlinear  function  of  the 
eigenvalue  Qx2  through  the  longitudnal  bending  and  twisting 
rigidity  terms  in  the  equation  for  the  basilar  membrane  as  well 
as  through  the  coupling  with  the  fluid.  It  is  also  complex 
valued  because  of  the  viscous  contribution  in  the  kinematic 
boundary  condition. 

For  use  in  the  possible  eigenfunction  expansions  to  come 
from  this  problem,  as  well  as  for  numerical  approximations,  it 


i 
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is  of  interest  to  know  the  integral  theorems.  The  simplest  is 
obtained  by  integrating  (4)  over  the  cross-section,  from  which 
one  finds  that 

_  _  G+ 

©x2(J7po  <3ydz  -  0/  po)  =  <*/  ?o  dy 

f  a*  g_ 


More  important  is  the  Rayleigh  quotient,  which  is  obtained  by 
multiplying  by  pg  and  then  integrating  to  obtain 


©x2</JP02  dydz  "  S/  Po2) 


3<r 


G+ 

G_ 


-  "  //[(PO  )2+(P0  ) 2  ]dydz  (13) 

¥  y  z 


+  (D^0x4-w2h  )  CQ2  Jdy 


This  is,  in  effect,  a  relationship  involving  the  kinematic  and 
potential  energies  of  the  system  with  the  viscous  dissipation. 
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SMALL  B  APPROXIMATION 


Even  though  boundary  layer  theory  was  used  to  obtain  (2) 
and  (8),  the  fact  that  B  is  small  was  not  considered  in  the  two- 
variable  expansions  discussed  in  the  last  section.  If  we  impose 
this  on  our  eigenvalue  problem  there  is  a  small,  but  important, 
simplification  in  the  boundary  conditions.  To  do  this  one  finds 
that  for  small  8 


ex  -  kg(x)  +  Bki(x)  + 


PO  ’  PQlX'Y'Z)  +  &Pl  +  ...  » 


c o  -  no(x,y)  +  Bni  +  . . .  , 


(14a) 


(14b) 


(14c) 


Introducing  these  into  (4)  the  following  eigenvalue  problem  is 


obtained 


(3y2+3z2)p0  =  k02p0  , 


(15a) 


9y2(D23y2ri0)  “  k02[ 3y2( vlD2no)  +  vlD23y2,iO  +  4 3y ( 3y  no >  J 


+  (ko4Dj.  -  u>2h)no  =  -2u>2po(x,y,0+)  , 


(15b) 


where 


^PO  = 


|  -a no  on  BM 
|  0  on  rigid  wall 


(15c) 


The  viscous  correction  kjjx)  in  (14a)  is  found  from  the  solva¬ 
bility  condition  for  the  0(B)  problem,  which  for  the  pressure  is 

( 3y2+3z2-kQ2)Pi  *  2kgkiPo 
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Multiplying  this  by  pg/  integrating  over  the  cross-section  and 
then  using  Green's  theorem  one  £inds  that 

*o/  PO2 

k.(x)=  _ £1 _  .  (16) 

1  G+ 

2//P()2dydz  ♦  Ig_  [D3(„0y)2*k0\„02}6y 

As  in  (7),  it  has  been  assumed  in  deriving  (16)  that  the  bending 
and  twisting  rigidities  are  only  functions  of  x.  It  is  also  of 
interest  to  note  that  this  result  can  be  obtained  from  the 
Rayleigh  quotient  by  substituting  (14)  into  (13).  In  any  case, 
the  complete  first  order  approximation  of  the  solution  is  then 
found  by  satisfying  (7)  along  with  the  boundary  conditions  (6), 
(11),  and  (12). 


SOLUTION  OF  NONLINEAR  EIGENVALUE  PROBLEM 

Now  that  the  original  system  of  partial  differential 
equations  has  been  reduced  using  slender  body  and  boundary  layer 
theory  it  remains  to  solve  the  nonlinear  eigenvalue  problem 
(15).  Once  this  is  done  it  is  a  simple  matter  to  find  k^  and 
satisfy  the  associated  boundary  conditions.  To  solve  (15)  one 
can  use  numerical  methods ,  modal  expansions  ,  or  Green's 
functions.  Each  has  its  advantages  as  well  as  its  restrictions. 
We  will  illustrate  the  solution  of  the  problem  using  the  last 
two  methods.  The  approach  is  similar  for  both  methods  since, 
for  a  given  geometry,  the  eigenvalue  problem  is  reduced  to  an 
integro-dif ferential  equation  for  nolx,y)  and  kg(x).  After  this, 
Fourier  expansions  are  used  to  reduce  the  problem  to  a  set  of 
algebraic  equations  that  have  kg2  as  the  eigenvalue  and  the 
Fourier  coefficients  as  the  eigenfunction.  To  facilitate  the 
presentation  we  now  assume  that  the  plate  is  homogenous,  so 
D2  “  h  =  1,  and  Dj_,  D3  are  constants.  However,  both  the  width 
of  the  basilar  membrane  and  the  transverse  cross  section  of  the 
cochlea  are  still  arbitrary  functions  of  x. 

Modal  Expansions 

For  this  case  we  assume  that  the  cross-section  is 
rectangular  (Fig.  2).  By  separating  variables  in  (15)  one  finds 
that 

m  ^ 

Po(x,y,z)  -  l  PmCoshxm(H-z ) •  cos-r^  ,  (17) 

m*0 


where 


-H(x) 


m — 

GJx) 


/////// iii _ 

H(x) 


FIG.  2.  Transverse  cross-section  used  for  modal 

expansion  solution  of  eigenvalue  problem. 


and 


for 


G+ 

>  =  ~acm  J  nocosymsds 

m  XmsinhXmH  J G_ 


m  *  0 

m  *  0 


Substituting  this  into  (15b)  one  obtains  the  following  integro 
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differential  equation  for  kg2  and  n< 


G+ 

( 3y4-2D3ko23y2+Dik04-a)2)  no  =  /  L(x,y;s)nods  , 

G_ 


where 


:,y;s)  =  I  ara(: 


L(x, y;s)  =  L  a  (x)cosYmy.cosY  s  • 
m=0  m  m  m 


(18) 


and 


a  =  2au)2cm 
m  Xmtanh  XmH 


To  solve  this  we  expand  no  in  modes  as  follows 

V*'*’  -  nLbnCOS<‘nly-(^)j  ' 


where 


(19) 


a  s  11  ( 2n-l )  m 
n  G+-G_ 

This  satisfies  the  simply  supported  boundary  conditions  (6),  and 
for  it  to  satisfy  (18)  one  finds  that  the  bn’s  and  kg2  have  to 
be  such  that 


(a  4+2D,k  2a  2+d  k  4_u2)b  =  _ 2 _  £  J  a  b  K  K  , 

4  30410  4  G+-G_  mZ0  m  n  mn  mt 


where 


G+ 

Kmn  -  iG_ooa’'m1'-cos“ 


]dy 


We  can  write  this  set  of  equations  in  matrix  form  as 


Mb  -  0  , 


(20) 


where  M  is  a  nonlinear  function  of  the  eigenvalue  kg2. 


i 


A 
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To  find  the  solution  it  is  necessary  to  know  how  many  modes 
are  needed  for  the  pressure  (17)  and  basilar  membrane  (19)  to 
obtain  a  reasonable  approximation.  If  the  basilar  membrane  is 
in  the  center  of  the  partition  so  that 

G+(x)  =  -G_ ( x )  =  G( x )  ,  (21) 

then  numerical  studies  show  that  the  first  mode  for  the  basilar 
membrane  and  the  first  five  modes  for  the  pressure  are 
sufficient.  However,  if  the  basilar  membrane  is  not  centered, 
or  if  the  problem  is  asymmetric  in  y  for  any  reason  (e.g., 
because  of  asymmetric  boundary  conditions)  then  it  is  necessary 
to  use  several  modes  for  the  basilar  membrane. 

For  the  case  of  a  centered  basilar  membrane,  the 
characteristic  equation  obtained  with  a  one  mode  approximation 
in  (19)  is 

<“r2t2D3"rk02tDlk04-“2)G  *  amKml2  >  1221 

m=u 

where 

“r  '  '  (23> 


and 


2  /  u)j»COS  y  mG 


if  Y  ^  *  w 
'm  r 


if  Ym2  =  «r 


This  problem  is  similar  to  the  one  obtained  by  Steele  and  Taber 
(1979)  using  the  method  of  averaged  Lagrangians. 


. 

■  WanW  ■  '  . 


A 
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Green's  Function 

For  this  approach  we  will  take  a  semi-circular 
cross-section  and  a  centered  basilar  membrane  (Fig.  3).  The 
Green's  function  G(x,y,z;y ' ,z ' )  satisfies 


(3v.2+3z.2-k02)G  =  6(y'-y)6(z'-z)  , 


where 


3 n i  G  -  0  along  S'?' 

•  •  fji 


(24) 


From  this  and  (15a)  it  follows  that 


G 

pQ(x,y,z)  =  2a/  G(x,y,z;y* ,0)nQ(x,y' )dy'  .  (25) 

-G 

To  find  G  recall  that  the  fundamental  solution  of  (15a)  is 
the  modified  Bessel  function  Kg,  and  so, 


G( x,y,z;y ' ,z ' )  =  ziKQ{ |kQ ^ (y-y • )2+(z-z ' )2)  +  g ( x,y , z ;y • , z ' )  , 

where  the  function  g  satisfies  (15a)  and  is  such  that  G 
satisfies  (24),  Therefore, 

„,v  v  =  K0'(k0R)I0(kor)Io(kor,> 

q(»,y,2.Y  - 

+  2  l  Kn'(k0R)In(lcor)In(lc0r’)  co3(ne) 
n=l  in ’  (*0^ 

where  r  *  /y2+z2  ,  r’  «  ^ (y ' ) 2+(z  ' )  ^  ’  and  ®  the  an9le 

between  (y',z’)  and  (y,z)  with  respect  to  the  origin. 


FIG.  3.  Transverse  cross-section  used  for  Green's 
function  solution  of  eigenvalue  problem. 

With  this  it  therefore  follows  that  no  satisfies  the  integro- 
differential  equation 

G 

(  3y4-2D3ko^3y^+Diko^-w^)  no  *  -2au)2  J  <( x , y ; s  )  no^s  •  (  26 

where 

K ( x,y; s  )  =  -  _i  KQ( |kQ  |  |y-s | )  +  g(x,y,0;s,0)  • 

It  is  possible  to  simplify  this  by  using  the  approximation 


kQR>>l,  which  is  the  case  of  interest  in  the  examples  to  be 
discussed  later.  Using  the  usual  approximations  for  Bessel 
functions  of  large  argument  and  the  identity 


IQ(w)  =  l  (-DnI  (S)In(o)ein*  ,  where  w  =  (S2+o2-2So  cos  <*>) 1/ 2 - 

—  oo 

the  eigenvalue  problem  becomes 


(8v4-2D3l'02sy2+Dlk04-“2ln0  ‘  ^ 


J  f (y-s ) nQds  ,  ( 27 ] 


where 


f  (s ) 


-2  Ik.  | R 


K0(  I^qs  |  )  +  tte  0  I0(  |k0s  |  ) 


(28) 


We  have  also  assumed  that  the  plate  is  centered  (21). 

For  the  first  mode  approximation  of  no/  one  finds  that  the 
characteristic  equation  for  kg^  is 

<Dr2+2D3U)rkg  2+D^ko^-ij2  =  U^G(  x  ) F  (  <q,  )  ,  (29) 


where 

1 

F(V<i>  =  —  /  [Kq(  (cQs)+Tre  <oKll0(  <qS)  J  [  ( 1-s  )cos  us+Asin  us  Jds 

(30) 

for  <g(x)  =  2|kg|G(x)  and  *i(x)  =  R(x)/G(x).  It  remains  to 
solve  this  equation  for  kg,  but  it  is  of  interest  to  note  that 
the  effect  of  the  outer  boundary  is  contained  entirely  in  the 
exponential  term  in  (30)  through  the  ratio  of  the  half  width 
of  the  basilar  membrane  to  the  cross-sectional  radius. 

L 


NUMERICAL  RESULTS 


To  illustrate  the  qualitative  nature  of  the  solution  let 
B  =  0.05  cm,  L  =  3.5  cm,  v  =  0.008  cm/sec,  and  p  =  1.0  gm/cm2. 
Also,  w=hcpc  and 


3 

D  *  =  E2hc 
2c  12(l-o2) 
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where  E2  =  4X10^  dyn/cm2,  hc  =  3.8  X  10-2  cm,  and  0  =  1/2.  The 
boundary  of  the  basilar  membrane  is  given  as 


G(x) 


__L( 5x+l ) 
12 


0<x<l  , 


and  the  cross-sectional  area  is  assumed  to  have  a  constant  value 
of  0.01  cm2,  which  determines  H  and  R  in  Figs.  2  and  3.  Although 
we  have  not  attempted  to  be  absolutely  realistic,  these  values 
are  representative  of  those  found  for  the  human  cochlea  (Holmes, 
1980). 

Before  discussing  the  full  solution  we  first  consider  the 
effect  of  the  bending  and  twisting  rigidities  of  che  plate,  as 
well  as  the  shape  of  the  cross-section,  on  the  eigenvalue  kg2. 

In  Fig.  4  the  values  obtained  using  (22)  are  shown  for  the  case 
of  an  isotropic  plate,  where  =  D3  =  1,  and  a  strongly  ortho¬ 
tropic  plate,  with  D3  =  D3  =  0,  for  a  driving  frequency  of 
1800  Hz.  Also  shown  are  the  values  obtained  for  a 
semi-circular  cross-section  for  the  strongly  orthotropic  plate. 

In  this  figure,  as  well  as  in  all  that  follow,  eight  pressure 
modes  are  used  in  (17)  and  (22).  From  this  result  one  can  see 
that  the  effect  of  an  orthotropic  plate  is  not  significant  near 


IJKil 


IKII 


IJtfll 


ir:ii 
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X-flXIS 


The  eigenvalue  kg(x)  as  a  function  of  longitudinal 
position  for  a  driving  frequency  of  1800  Hz.  Shown 
are  the  values  for  a  rectangular  cross-section  (22) 
with  03=03=1,  lower  curve,  and  with  03*03=0,  upper 
curve.  Also,  shown  is  the  solution  for  a  circular 
cross-section  (29)  with  03=03=0,  which  is  essentially 
indistinguishable  from  the  analagous  curve  for  a 
rectangular  cross-section. 
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the  basal  end,  x  =  0,  as  the  wavenumber  is  still  relatively 
small.  However,  as  the  wavelength  decreases,  so  x  increases, 
the  longitudinal  bending  terms  have  a  substantial  effect  on  the 
eigenvalue.  In  other  words,  the  smaller  the  wavelength  the  more 
important  the  longitudinal  coupling  is  in  the  basilar  membrane. 
This  difference  also  changes  the  damping  in  the  sense  that  the 
smaller  the  eigenvalue  the  less  the  viscous  dissipation  of  the 
wave  (16).  As  for  a  rectangular  and  circular  cross-section, 
the  curves  in  the  figure  are  indistinguishable.  Equation  (29) 
was  used  to  calculate  the  semi-circular  geometry,  but  only  for 
k0>l,  as  was  assumed  in  the  derivation. 

There  have  been  a  number  of  studies  on  the  elastic 
properties  of  the  basilar  membrane  which,  unfortunately,  usually 
involve  relatively  large  amplitude  deformations  (von  Bekesy, 
1960;  Voldrich,  1978).  However,  it  appears  that  the  longitud¬ 
inal  coupling  is  relatively  weak,  and  so,  in  the  next  examples 
it  is  assumed  that  the  basilar  membrane  is  strongly  orthotropic, 
which  means  that  =  D3  =  0.  In  this  case  it  is  a  simple 
matter  to  prove  that  kg2  is  real  since,  from  (13),  one  finds 
that 

9  G  + 

k02  / /  I P 0  I  dVdz  =  “/ /  I VP0 1 2  dydz  +  “  /  |n0  |2  dy  . 

u  »  u  f  u  2  o)2  G-  uyy 

In  the  examples  discussed  here  there  is  only  one  positive 

eigenvalue,  which  leads  to  left  and  right  traveling  waves  in 

(3).  When  the  viscous  correction  is  included,  i.e.,  kj_, 

these  waves  are  damped.  For  the  parameter  values  given 

above,  when  the  frequency  is  above  about  600  Hz,  it  is  not 
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necessary  to  include  the  left  running  wave  as  the  apical 
boundary  condition  is  effectively  satisfied  by  the  one 
traveling  to  the  right.  In  the  examples  presented  below  the 
reflected  wave,  that  is,  the  one  moving  towards  the  stapes, 
is  not  included. 

With  the  assumption  that  the  reflected  wave  is  negligible 
the  solution,  as  given  by  (3),  is  relatively  easy  to  find.  To 
do  so  it  is  also  assumed  that  the  problem  is  symmetric  about  the 
z-axis,  so  that  it  is  sufficient  to  only  use  the  first  mode  in 
(19).  In  this  case  the  displacement  of  the  basilar  membrane  is 

i  [t-*/2-k(x,  e  )  ] 

n*(x,y,t)  -  As * •  A ( x ) e  cos(iry/2G)  ,  (31) 


where 


k(x,e)  «  I  /  [kQ(s)+6Qk1(s) ]ds  , 

£  0  u  u  x 


An  -  L  k, (S )dS 

A  ( x )  =  _ _ A0  ,  •  e  *  0  0  1 

McToO/jpT^dydz 

¥ 


Ao  -  H  P010.y,z)2dydz 


0O  =  i/n  ,  Aw  =  AwVbZ  is  the  (dimensionless)  area  of  the 
oval  window,  and  As*  is  the  amplitude  of  the  stapes.  The  func¬ 
tion  kfl(x)  is  determined  from  (22)  for  a  rectangular  cross- 
sectional  geometry  and  from  (29)  for  a  semi-circular  cross- 
section.  After  this,  Po(x>y/z)  is  found  from  either  (17)  or 
(25)  and  ki(x)  is  determined  from  (16).  It  is  of  interest  to 
note  that  (31)  is  equivalent,  in  terms  of  its  functional  form,  to 
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the  solutions  obtained  by  Steele  (1980)  and  by  Lighthill  (1983). 
However,  their  procedure  to  calculate  the  coefficients  is, 
apparently,  different  so  it  is  unclear  how  their  solutions 
compare  with  ours. 

The  phase  ph(x)  of  the  wave  in  (31)  is  given  as 

x 

ph  ( x )  =  -I  -  A  /0[kfl(s)  +  60k:1(s)  jds  .  (33) 


The  phase  velocity  v*(x)  is  therefore 


v*(x)  =  _ _ ^ _ 

k0(x)+eo*i(x ) 


(34) 


The  function  A(x)  in  (32)  is  the  amplitude  of  the  wave  on 
the  basilar  membrane.  To  see  how  it  depends  on  x  note,  from  Fig. 
4,  that  kg  increases  significantly  with  x.  For  large  kg  one  can 
approximate  the  eigenvalue  problem  (29)  with 

D1k()4  +  2D3<ork02  +  u,r2  -  u>2  =  2au2  .  (35) 


Also,  for  large  kg  in  (16) 


*1  -  kQ2 


(36) 


and,  from  (17), 


(37) 


If  Di»D3*0  then 
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k 


0 


2ati)2 

aj  2_u2 

r 


whereas,  if  0^=032^0  then 


k 


0 


(38) 


(39) 


This  means  that  for  the  strongly  orthotropic  case  (38),  kg 

2 

becomes  infinite  as  x  increases,  that  is,  as  w  =  (__! ) 

r  2g' 

approaches  u>.  This  represents  the  position  of  the  transverse 
beam,  in  the  basilar  membrane,  with  a  fundamental  frequency 
corresponding  to  the  given  driving  frequency.  For  an  ortho¬ 
tropic  case  as  in  (39),  kg  has  no  singularity  and  simply 
increases  continuously  with  x.  The  singular  nature  of  the 
strongly  orthotropic  case  results  in  a  nonuniform  expansion  in 
(14).  However,  from  (32)  and  (36),  the  solution  decays  exponen¬ 
tially  as  it  enters  this  region.  In  fact,  the  place  principal 
arises  from  the  balance  in  (32)  between  the  increase  in  the 
coefficient,  as  given  in  (37),  and  the  exponential  decay  coming 
from  (36)  and  (38).  As  the  wave  propagates  down  the  basilar 
membrane  the  phase  velocity  decreases,  and  goes  to  zero  as  it 
approaches  the  singular  point.  This  is  reminiscent  of  the 
behavior  associated  with  critical  layer  absorption  as  discussed 
by  Lighthill  (1981).  It  is  also  of  interest  to  note  that  the 
large  wavenumbers,  as  determined  from  (35),  are  altogether 
independent  of  the  shape  of  the  cross  section.  At  the  same  time, 
the  resonant  point  (38)  would  not  exist  without  the  fluid  or 
without  the  mass  of  the  basilar  membrane. 


ft 
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In  Fig.  5  the  amplitude  and  phase  of  the  wave  propagating  down 
the  center  of  the  basilar  membrane  (y=0)  is  shown  for  three  driving 
frequencies  (Appendix  C).  One  can  see  that  the  amplitude  at  600  Hz 
is  not  negligible  at  x  =  1,  which  shows  that  it  is  necessary  to 
include  the  reflected  wave  for  lower  frequencies.  Also,  the  ampli¬ 
tude  for  all  three  frequencies  is  not  zero  at  x  =  0,  which  is  the 
result  of  not  including  the  contribution  of  the  edge  layer.  If 
this  were  done  there  would  only  be  a  minor  change  in  the 
curves,  except  for  the  very  high  frequencies  which  have  a 
maximum  near  to  this  region.  In  any  case,  the  amplitude  curves 
show  the  usual  shift  towards  the  basal  end  as  the  frequency 
increases.  Also,  the  amplitude  increases  from  75  at  600  Hz  to 
220  at  6000  Hz.  In  other  words,  the  maximum  amplitude  of  the 
basilar  membrane  is  from  115  to  225  times  the  amplitude  of  the 
stapes  for  frequencies  between  1000  and  10,000  Hz,  which  is  con¬ 
sistent  with  the  measurements  of  Rhode  (1971)  and  Khanna  and 
Leonard  (1981).  It  is  also  interesting  to  note  that  _in  situ  the 
amplitude  of  the  stapes  decreases,  for  a  given  sound  pressure, 
as  the  frequency  increases,  which  would  tend  to  counteract  this 
growth . 

Tuning  and  phase  curves  are  shown  in  Fig.  6,  which  represent 
the  amplitude  of  the  wave  at  a  fixed  location  as  a  function  of 
frequency.  The  points  x  =  0.82,  0.45,  0.20  are  chosen  as  they 
represent  the  positions  of  the  maximum  amplitudes  in  Fig.  5. 
Overall,  the  curves  show  a  relatively  sharp  tuning.  Whether  or 
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FIG.  6.  The  amplitude  and  phase  of  the  wave  on  the  basilar 

membrane,  as  a  function  of  frequency.  The  three  fixed 
spatial  locations  shown  correspond  to  the  points  of 
maximum  amplitude  in  Fig.  5. 
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not  they  are  physiologically  realistic  is  unclear  but  they  are 
significantly  sharper  than  those  obtained  by  von  Bekesy  (1960). 
However,  they  are  not  as  sharp  as  would  be  implied  from  the 
neural  tuning  curves. 

The  phase  velocities  of  the  waves  in  Figs.  5  and  6  are 
shown  in  Fig.  7,  along  with  von  Bekesy 's  (1960)  measurements. 

In  Fig.  7a  it  is  seen  that  the  wave  slows  down  as  it  propagates 
down  the  basilar  membrane,  and  significantly  so  once  it  begins 
to  approach  the  point  of  maximum  amplitude.  The  values  found 
from  the  theory  are  smaller  than  those  obtained  by  von  Bekesy, 
which  he  states  (pg  458)  "are  representative  of  all  the  measure¬ 
ments"  made  on  cochlear  models  and  cadevar  specimens.  From  Fig. 
7b,  at  a  fixed  spatial  location,  the  phase  velocity  has  a  rela¬ 
tive  constant  value  until  the  frequency  begins  to  approach  the 
frequency  of  largest  amplitude  for  this  point.  When  this  occurs 
the  phase  velocity  decreases  very  rapidly  as  the  frequency 
increases . 

It  is  rather  easy  to  determine  the  position  x^  of  the  max¬ 
imum  amplitude  as  a  function  of  frequency.  The  result  is  shown 
in  Fig.  8,  which  also  includes  the  measurements  of  Crowe,  et  al. 
(1934)  on  the  spatial  localization  based  on  pathological 
condition,  as  well  as  von  Bekesy 's  (1960)  measurements  on  cada¬ 
vers.  It  appears  that  in  this  aspect  the  theory  agrees  well 
with  experiment. 

The  last  result  concerns  the  difference  limen  (DL),  or 
just-noticeable-differences,  obtained  from  the  theory.  These 
represent  the  minimum  change  in  the  driving  frequency  for  there 


VELOCITY  (CM/SEC)  VELOCITY  (Ctt/SEC) 


FREQUENCY  (HZ) 
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FIG.  8.  The  longitudinal  location  x^  of  the  maximum  amplitude 
as  a  function  of  frequency  as  determined  from  the 
theory  (solid  curve).  Also  shown,  by  the  horizontal 
bars,  are  the  spatial  locations  measured  by  Crowe,  et 
al  (1934),  and  von  Bekesy's  (1960)  measurements  are 
given  by  the  x's. 
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to  be  a  perceivable  difference  to  an  observer.  It  should  be 
kept  in  mind  in  the  discussion  that  such  psychoacoustic 
measurements  are  subject  to  rather  significant  individual 
differences  as  well  as  influenced  by  such  things  as  resonances 
of  the  ear  canal  (Henning,  1966).  In  accordance  with  the  place 
principle  of  hearing,  we  assume  that  the  perceived  change  in 
frequency  is  due  to  the  shift  of  the  maximum  amplitude  of  the 
wave  on  the  basilar  membrane.  The  question  is  how  far  does  it 
have  to  shift?  The  answer  is  related  to  the  stimulation  of  the 
hair  cells.  There  are  approximately  3500  rows  of  such  cells 
along  the  membrane,  and  we  assume  that  the  maximum  point  has  to 
shift  at  least  the  distance  of  one  row  before  it  is  possible  to 
recognize  a  difference.  The  resulting  difference  limen  curve  is 
shown  in  Fig.  9  along  with  the  experimental  results  of  Wier,  et 
al  (1977),  Nordmark  (1968)  and  Shower  and  Biddulph  (1931). 

Also  included  is  the  curve  one  would  obtain  if  the  shift 
required  is  three  rows.  All-in-all  the  agreement  between  theory 
and  experiment  is  very  good.  It  implies  that  the  system  is 
responsive  to  very  small  shifts  in  the  frequency,  and  the 
spatial  difference  is  approximately  constant  between  1000  and 
10,000  Hz.  The  theoretical  curves  in  the  figure  can  be  written 
as 

log  (DL)  =  0. 023 /ft  +  bo  ,  (40) 

where  bQ  =  -1.08,  -0.778,  -0.602,  -0.477  for  shifts  of  one, 
two, three,  four  hair  cell  rows,  respectively.  Interestingly, 
the  linear  relation  (40)  has  been  obtained  from  empirical 
curve  fits  by  Wier,  et  al  and  also  by  Nelson,  et  al  (1983). 
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.  The  difference  limen  (DL)  obtained  from  the  theory  for 
a  spatial  shift  of  the  maximum  amplitude  corresponding 
to  a  single  row  of  hair  cells  (lower  solid  curve)  and 
to  three  rows  of  hair  cells  (upper  solid  curve).  Also 
shown  are  the  experimental  measurements  of  Shower  and 
Biddulph  (1931),  with  the  +,  Nordmark  (1968),  with  the 
x,  and  Wier,  et  al  (1977),  with  the  0. 
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The  latter  show  that  the  slope  in  (40)  is  relatively  constant  as 
the  stimulus  sensation  level  (SL)  increases,  whereas  the  intercept 
bQ  decreases.  For  SL's  between  10  and  80  dB  the  slope  for  the 
data  of  Nelson,  et  al,  is  0.0214,  for  Wier,  et  al,  it  is  0.0238, 
and  for  Harris  (1952)  it  is  0.0235.  Clearly,  all  three  slopes 
are  very  close  to  our  value  of  0.023.  As  for  the  decrease  of  bQ 
as  the  SL  increases,  this  would  imply  from  the  theory  that  a 
smaller  spatial  shift  is  necessary  to  discern  a  difference  in 
frequency  for  a  higher  SL.  This  is  consistent  with  the  fact 
that  the  amplitude  on  the  basilar  membrane  also  increases 
with  the  SL,  which,  in  turn,  means  an  increase  in  the  neural 
signal  from  the  hair  cell. 


k 
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SUMMARY 


Using  slender  body  and  boundary  layer  theory  we  have 
reduced  the  three-dimensional  hydroelastic  model  for  the  cochlea 
to  a  nonlinear  eigenvalue  problem  in  the  transverse 
cross-section.  This  problem  comes  from  the  strong  coupling 
between  the  basilar  membrane  and  fluid,  and  does  not  exclude  any 
of  the  effects  of  the  plate,  such  as  the  inertia.  The  subse¬ 
quent  response  of  the  basilar  membrane  to  a  pure  tone  of 
sufficiently  high  frequency  is  a  right  running  wave  that  begins 
to  grow  in  amplitude.  However,  the  wavelength  decreases  and  so, 
due  to  the  viscosity  of  the  fluid,  the  wave  eventually  reaches 
a  spatial  location  after  which  it  quickly  decays  to  zero.  For 
lower  frequencies  the  reflected  wave  is  not  negligible 
particularly  in  the  region  near  the  helicotrema,  and  has  to  be 
included  in  the  solution.  This  wave  would  be  even  more 
important  if  there  was  no  damping,  as  the  response  to  a  pure 
tone  would  be  simply  a  standing  wave  and  one  would  obtain 
resonances  corresponding  to  the  eigenvalues  of  the  system 
(Holmes,  1979). 

As  outlined  at  the  beginning,  there  are  a  number  of 
assumptions  made  to  formulate  the  model.  However,  from  the 
resulting  analysis,  it  has  been  shown  to  describe  a  good  deal  of 
the  qualitative  behavior  known  to  occur  in  the  cochlea. 
Nevertheless,  there  are  a  number  of  interesting  extensions  of 
the  analysis  which  could  be  undertaken.  For  example,  the 
maximum  amplitude  on  the  basilar  membrane  at,  say,  1800  Hz  is 
about  160  times  larger  than  the  amplitude  of  the  stapes.  This 
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raises  the  question  as  to  whether  it  is  necessary  to  consider  a 
nonlinear  hydroelastic  theory  in  the  neighborhood  of  this  point 
to  account  for  these  relatively  large  amplitudes.  Also,  there 
has  been  a  question  as  to  where  the  damping  in  the  system  is. 

We  have  only  considered  the  dissipation  due  to  the  viscosity 
of  the  fluid,  but  there  surely  is  damping  in  the  basilar 
membrane  and  its  associated  structures.  This  is  particularly 
true  for  the  short  waves  which  are  present  on  the  basilar 
membrane  on  the  apical  side  of  the  maximum  amplitude.  It  is 
not  clear  exactly  how  this  damping  manifests  itself  but  most 
of  the  standard  viscoelastic  theories  can  be  included  in  our 
analysis  without  increasingly  the  difficulty  significantly. 


V-  *>r«V  .  ■  _ 
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APPENDIX  A 

Derivation  of  Boundary  Condition  for  the  Pressure 

The  perilymphatic  fluid  in  the  upper  chamber  is  assumed  to 
be  Newtonian.  The  equations  of  motion  come  from  the  linearized 
Navier-Stokes  equations  and  the  incompressibility  condition 
which,  in  scaled  coordinates,  are 

( i-62y2)v  =  -7p  ,  (Ala) 

V*v  =  0  .  (Alb) 

To  facilitate  the  presentation  we  have  let  e  *  1  in  (1).  Also, 
the  time  periodic  solution  is  being  considered  here. 

It  follows  immediately  from  (Al)  that  p  satisfies  (la). 

As  for  the  boundary  condition,  it  is  assumed  that  there  is  a 
prescribed  normal  velocity  component  as  in  (2).  To  give  it 

4 

explicitly  let  xq  be  a  point  on  the  boundary  surface  S,  which  is 

4 

assumed  to  be  regular  in  a  neighborhood  of  xq-  It  is  then 

4 

possible  to  introduce  local  cartesian  coordinates  at  xo  of  the 
form 

4-  4  4  4  4- 

x  -  x0  =  +  s2e2  +  nen  , 

♦  +  + 
where  en  represents  the  unit  outward  normal  to  S  at  xq,  and  ej, 

♦  4  4  4  4 

e2  are  in  the  tangent  plane  to  S  at  xq.  Note  en,  ei,  e2  are 

4 

independent  of  the  three  spatial  coordinates.  The  velocity  v 
can  now  be  written  as 

v  -  vs  +  vnen 


JaiLL 
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where  vn  and  vs  are  the  normal  and  tangential  components  of  v, 
respectively.  The  gradient  operator  can  be  similarly  factored 
as 

7  =  7  +  e  _i  • 
s  n  3n 


With  this,  the  boundary  condition  at  xg  is 

v(xg)  =  -inen  .  (A2) 

The  fluid  motion  in  the  inviscid  region  is  found  from  (Al) 
by  letting  6+0  with  x  fixed  in  the  fluid.  The  appropriate 
expansions  in  this  case  are 

v  "  vo(s^,S2»n)  +  6vi  +  ...  , 

and 


As  is  usual,  it  is  impossible  to  satisfy  (A2)  with  this 
inviscid  approximation,  and  so,  it  is  necessary  to  consider  the 
viscous  boundary  layer  on  S.  With  the  boundary  layer  coordinate 
n  *  n/6,  the  expansions  for  this  region  are 

V  ~  Vg(s^,S2»n)  +  6V^  +  ...  ,  (A3) 


Introducing  these  into  (Al)  one  finds  that  the  first  terms  in 
these  expansions  satisfy 


Solving  these  and  matching  the  solutions  with  the  first  terms 
for  the  inviscid  region 

?0  =  P0<sl's2'°>  ' 

Vn0  *  -in(si,s2)  r 

and 

Vs0  *  vs0(slfs2,0)(l-e/i  n)  . 

In  matching  these  solutions  it  is  assumed  that  the  fluid  lies  in 
the  region  n<0. 

Now,  from  (Al ) , ( A3 ) , { A4 ) 
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So,  solving  and  matching  with  the  outer  solution 
P1  =  n3np0(si,s2»0)  +  p1(s1,s2,0) 


and 

vnl  =  ~7s*vsO(s1'S2,0)  [n+d-e-1^  n)//T] 

=  3nvn0(sl's2»°> [n+( 1-e"1^  n)//T]. 

From  this,  for  the  inviscid  expansion  to  match  with  the  one  for 
the  boundary  layer,  vn  must  be  such  that 

vn(si,s2,0)  -  -in  +  63nvno(si,s2,0)  +  ... 

Putting  our  results  together  then 

P(si,s2,n)  '  Po(si,s2,0) 

+  <5  [  n3nPo(si,s2,0)  +  pi(si,s2,0)]  +  ...  (A5) 


and 


3nP(sl's2'°>  =  -ivn(si,s2,0) 

=  -n  +  B  3n^Po (sl ,s2 , 0 )  +  0( 62) 
=  -n  +  03n2p(s1,s2,O)  +  0(62). 


Therefore , 


(3n-B3n2)P  *  -h  on  S. 


Also,  from  (A5 )  we  obtain  the  pressure  on  the  basilar  membrane 
used  in  (16). 
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As  is  clear  from  the  derivation,  these  expressions  are 
valid  up  to  0(62).  Moreover,  we  have  used  a  tangent  plane 
approximation  of  the  surface  near  xg«  Given  the  local  nature  of 
the  boundary  layer  this  is  sufficient  but,  if  a  global  analysis 
is  desired,  one  can  use  orthogonal  curvilinear  coordinates  as 
presented  by  Howarth  (1959). 
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APPENDIX  B 

Solvability  Condition 

To  derive  the  solvability  condition  (7),  which  gives  the 
slow  modulation  of  the  wave,  we  assume  here  that  Dj_,  D2,  are 
constants.  In  this  case  the  0(e)  problem  obtained  from 
substituting  (3)  into  (1)  is 

Oy2+az2-0x2)pi  =  -i  (©xxP0  +  2ex3xP0)  /  (Bl) 

P<C1>  ~  iBl  (  40  x^  3  xc  0+^®x2®xx  5  0 )  2iD3  3y2 ( 20x  3X  S0+Qxx  £(}) 

=  -2u>2pi{x,y,0)  ,  (B2 ) 


where  P(  • )  represents  the  differential  operator  applied  to  cq 
in  the  0(1)  problem  (4b).  The  boundary  condition  from  (2)  is 


OnT-»3nT2)Pl  +  iOx[bPO-@<2b+PO>3nTPoj 


(B3 ) 


and  the  simply  supported  boundary  conditions  for  the  plate  give, 
for  y  =  G+(x)  , 

tl  =  3y2ti  “  2i®x^xD3 3y  50  =  ®  •  (B4) 


Now  multiplying  (Bl)  by  pg  and  integrating  over  the 
cross-section 


//  Po(V2-ex2)Pidydz  =  -ijj  3X( ©xPo2 )dydz 
¥  ¥ 


'  _i3x  i!-  SxP 0 -  i6x  Lbp0 


3  *f 


l 


However,  using  Green's  theorem 


If  p0{v2-ex2)Pldydz  =  /  (Po3n_Pl  -  Pl3nJ?0> 
t  t  t 

_  _  G+  _  _ 

=  -iex  /  [bp02-63  (bp02)]  +  af  ( CiPO” ^oPl >dy  » 

d*  T  g_ 


where  we  have  also  used  (B3)  and  (5).  Combining  this  with  the 
previous  result 


G+ 

■  zi*  j 

2oj2  JG. 


Ul 


<<0)_C0 


Ul)+2iDl3x(0x 


—  2iD3 Cg dy2 ( 2©x 3X Cg+ ®xx ^0 ^  Jdy 
G+ 

■  "  J  af  U.  [Dl0x3tO2+D36x<8y‘o>2Jdy  ' 

where  (B4)  has  been  used  to  obtain  the  last  equality.  From  this 
one  immediately  obtains  (7). 
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APPENDIX  C 

Solution  for  a  Rectangular  Cross-Section 

The  displacement  of  a  simply  supported,  homogenous, 
orthotropic,  and  centered  basilar  membrane,  as  determined  by  the 
right  running  wave,  is 

n*( x,y , t )  ~  A  *A(x)e1Lt-k(x'e)_  "/2jCos  (IX)  ,  (Cl) 

s  or. 


where 


k(x,£)  =  A  /  [kgfsj+ggk^s)  ]ds 


( C  2 ) 


"  -  M  kl(s)ds 

A( x )  =  •  e  e  0  0  1  ,  ( C 3 ) 

/k0<x)ii  Pq^  dydz 


TFAwkQ(  0) 

4G(  0) 


JkQ(0) J  fpQ(0,y,z)2  dydz 


( C4  ) 


The  parameters  are  (Figs,  1,2) 


e 


f 


a 


PB 

V 


and 


u  =  2irf  L . 

^  D 

Here  f  is  the  driving  frequency  (in  Hz),  As*  is  the  amplitude  of 
the  stapes,  and  Aw  *  (Aw*/B2)  is  the  dimensionless  area  of  the 
oval  window. 
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For  the  case  of  a  rectangular  cross-section,  the  eigenvalue 
)Cq(x)  is  found  by  solving  the  equation 


where 


4  .  or^.1.2  ,  2  2  _  i  J  „  2 

ui  —  u)  —  l  a  k 


°lk0  +  2D3“rk0  +  T 


G  m=0  »  ml 


and 


m 


u> 

r-y>2  ' 

a  =  . 

2  au^Cjn 

m 

Amtanh 

-  / 1,  2  . 
K0 

+  v  2  ,  Y  = 

’m  'm 

m  h 

H 

9 

) 

( — 

}  2H 

if 

m=0 

1  i 

H 

V 

if 

m*0 

1 

(  2/u^T  coSYmG 

if 

2 

Kml  = 


u  -y  2 
r  'm 


m  r 


if  Ym2=wr 


In  these  expressions,  and  D2  are  constants  but  G  and  H  are 
functions  of  x  (as  determined  by  the  geometry). 

Once  kg  is  found,  the  function  kj(x)  in  (C2)  and  (C3)  is 
found  from  (16).  To  do  this  note,  from  (17),  that 

P0<x*y*2)  *  I  Pm  cosh Am ( H-z ) cos rmy 
m=0 


9 


(C5 ) 
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where 


P 


m 


aCm^ml 
Am  sinhAmH 


Thus , 


k1  (x)  =  kQ(x) 


•?  «  n+m_ 

H  Z  p  “■(  1+cosh"  A  H)  +  Z  Z  (-1)  ppL 
m=0  m  ra  m=0  n=0  m  n  mn 


1  -2 - -2 


H  1 
m  = 


1  p  2(H  +_L_  s inh 2H A  )  +  2a(D1kn2  +  D,w  )G 
=  0  m  2Am  m ;  0,2^  1  0  3  r  1 


where 


s inh ( An+ Am )H  +  s inh ( An- Am )H  ^  n 


*n+*m 


*n-  *m 


Jmn 


H  +  _±_  sinh2 A  H 


2  A 


n 


n 


if  n=m 


The  denominator  of  the  coefficient  of  the  exponential  in 
(C3 )  can  be  calculated  from  the  relation 


//  p  2  dydz  =  I  H  l  p  2(l  +  ,  3  sinh2A  H] 

This,  at  the  same  time,  determines  the  constant  Aq  in  (C4). 

The  displacement  of  the  basilar  membrane  is  completely 
determined  from  the  above  formulas.  As  for  the  fluid  pressure 

p*(x,y,z,t)  -  As*(4ir2uf2)A(x)ei  [t-k(x,  e)-n/2  ]p0(x,y,z) 

where  Po(x»y*z)  is  given  in  (C5).  From  this  one  can  calculate 

the  velocity  v*(x,y,z,t)  of  the  fluid,  away  from  the  immediate 
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vicinity  of  the  cochlear  wall  with  the  relation 


* 

v  = 


*  * 


2  TT  f 


7  p 
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